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model with correlation is considered. Explicit formulas for the probability maximizing 
function and the cost reduction function are derived. Applicability of the results for 
the widely traded derivatives as digital, quantos, outperformance and spread options 
is shown. 
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1 Introduction 

> 

As recent events on the market have shown the risk appearing in pricing of financial 
contracts should be more thoroughly surveyed. Although the problem of minimizing risk 
ly-j ■ is widely studied in the literature, the great majority of the results do not meet the 

expectations of practitioners who are interested in straightforward applications. This 
paper is concerned with the issue of risk analysis for the multiple assets derivatives and 
provides explicit computing methods for the risk parameters. 

The risk is measured by the possibility of a partial hedging of the payoff. Thus our 
approach is based on the idea of quantile hedging which was introduced in [6j and later 
developed in various directions, see for instance [3], [10], [2], pQ. Let us briefly sketch a 
general concept. Denote by H a contingent claim and assume that the arbitrage free pricing 
method indicates its price p{H). This means that if the investor has an initial endowment 
x > p{H) then he is able to follow some trading strategy such that his portfolio hedges H 
with probability 1. If this is the case, then x carries no risk and the probability maximizing 
function $i equals 1, i.e. <&i(x) = 1. On the other hand, if x < p(H) then the shortfall 
probability is strictly greater than zero for each trading strategy and then 3>i (x) < 1. The 
grater the probability of shortfall is the smaller the value is. Thus the function <i>i 

can be viewed as a measure of the risk sensitivity to the price reduction of the option. 
There is also another aspect of the problem. Assume that the hedger is willing to accept 
some risk measured by the shortfall probability in order to reduce initial cost. He chooses 
a number a S [0, 1] and searches for a minimal initial capital $2(0) which allows to find a 
strategy such that the probability of the shortfall is smaller then 1 — a. Thus if the hedger 
accepts no risk, i.e. a = 0, then the minimal cost required to replicate H is just p(H). 
In this case the cost reduction function satisfies $2(0) = p{H). However, if a > then 
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<3?2(o0 < p{H) and the function $2 enables us to view the effect how the risk acceptance 
affects the cost reduction of the option. Recall the numerical example from |6J p. 261 
which shows that $2(0,05) = 0,59 • p{H) for a call option with certain parameters. This 
means that the acceptance of a 5% margin of risk reduces the hedging cost by 41%. This 
shows that quantile hedging is an attractive tool for the risk analysis and should be taken 
into account by traders. 

The basic problem, however, is to determine functions and $2 for specific deriva- 
tives. There are only a few examples in the literature where they are explicitly found. 
In [6] explicit formulas are given for the most important case of a call option in a clas- 
sical Black-Scholes model. The method can be mimicked to obtain formulas for the put 
option. The idea is based on reducing the original dynamic problem to the static one 
which can be solved with methods used in the theory of statistical tests. Since the market 
was complete the solution of the static problem could be obtained, via Neyman-Pearson 
lemma, by indicating a non-randomized test for the appropriate probability measures. 
The Neyman-Pearson lemma can be generalized for the case of composite hipotheses, i.e. 
when measures are replaced by the families of measures, see [1] where the solution in the 
abstract form is presented. However, straightforward applicability of this result towards 
incomplete markets seems to be questionable. This paper is devoted to determining func- 
tions $1 and $2 for the multiple assets derivatives in the Black-Scholes framework with 
correlation. As the market is complete, we follow the same general method as in [6], but 
we find the solutions explicitelly using specific features of the model. More preciselly, we 
show that the original problem can be reduced to that of finding another two deterministic 
functions ^1,^2 depending on H, which turned out to be regular, i.e. continuous and 
strictly monotone if H is of a reasonable form, see Propsosition 13.41 and Proposition 13.51 
Then, roughly speaking, $1 = f 1 cj/j 1 and $2 = ^2°^i 1 ', for a precise formulation see 
Theorem 13.61 In the one dimensional case when H is a call option the result covers the 
above mentioned example from [6]. We also determine explicit forms of and ^2 for 
commonly traded derivatives, see Section 2] and its subsections. As ^2 are rather of a 
complicated form, the inverse functions can not be given by analytic formulas but can be 
determined with the use of numerical methods. Thus a great advantage of our results is 
that they can be used in practice. 

The paper is organized as follows. In Section [5] we briefly recall the multidimensional 
Black-Scholes model and formulate the problem strictly. Section [3] contains the main result 
- Theorem 13.61 which is proceeded by a general discussion on the results from [6] and the 
Neyman-Pearson technique. The method established in Theorem 13.61 is used in Section 0] 
for calculating the functions $1, ^2 for two assets derivatives which are widely traded. At 
the beginning of Section 0] we adopt general setting to a two dimensional model and recall 
basic properties of a multidimensional normal distribution. In the sequel explicit formulas 
for $1 and ^2 are presented in the case of digital option, quanto domestic, quanto foreign, 
outperformance and spread options. 
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2 The model 



Let (f2, Ttyt G [0, T], P) be a fixed probability space with filtration. The prices of d shares 
are given by the Black-Scholes equations 



dS\ = S l t {otidt + a i dW i 



tii 



l,2,...,d, t€ [0,T], 



where a* G R, <7j > 0, i = 1,2, and W 4 = (W/, , ...,W/),i G [0,T], is a sequence 
of standard Wiener processes adapted to {.F^i G [0, T]} with the correlation matrix Q of 
the form 



Q 



1 Pl,2 Pl,3 
P2,l 1 P2,3 

Pd.l Pd,2 Pd,3 



Pl,d 
P2,d 



where 



Pi,j 



cor 



{wi,w(} 



i,j = 1,2, ...,d. 



We assume that Q is positive definite. The proces given above will be called a Q- Wiener 
process. The trader can invest his money in stocks as well as put it on a savings account 
which dynamics is given by 

dB t = rBtdt, t G [0,T], 
with r standing for a constant short rate. 

Remark 2.1 The most common approach for the description of the market is based on 
a sequence of independent Wiener processes, see for instance a classical textbook J^j. It 
can be shown that the model described above is equivalent to the model with d independent 
Wiener processes and the dxd diffusion matrix with constant coefficients. We work with a 
correlated Wiener process because it is more convenient for later calculations. Let us also 
mention that parameters in such model can be easily estimated from data, see p.104- 

Let us now briefly characterize a martingale measure of the model, i.e. a measure P 
which is equivalent to P such that the discounted price processes S\ := e~ rt S\,i = 1, 2, d 
are local martingales. The following is a version of Theorem 10.14 in [3] adapted to our 
finite dimensional setting. 

Theorem 2.2 Let ip be a predictable process taking values in W 1 satisfying 
Then the process 

W t = W t - fq^sds, t G [0, T], 
Jo 

is a Q- Wiener process with respect to the measure P with a density 



dP 
dP 



e J- T (Q-? V udWt)-% !o\ft?dt_ 
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It can be shown that each measure equivalent to P can be characterized by a density 
process 



Zi 



t€ [0,21, 



(2.1) 



for some predictable M d - valued process cp. The process S l is a P local martingale if and 
only if S l Z is a P local martingale. Thus the measure P can be determined by finding a 
process (p in (|2.ip such that S l Z, i = 1,2, d are P local martingales. Simple calculations 
based on the Ito formula yield 



ft 



Q 



a - rl. 



a 



- cti—r 

a.2—r 
o~2 



te[0,T]. 



The martingale measure P is thus unique and given by the density process 

Z t := e-W-M^mKlQ-^^]^ t € [0 ,r]. (2.2) 
Moreover, it follows from Theorem 12.21 that the process 

W t := W t + a - rU t, t£[0,T], 
a 

is a Q- Wiener process under P. The dynamics of the prices under the measure P can be 
written as 

dS\ = S l t (rdt + (XidWl) , i = 1, 2, d. 
The wealth process with the initial endowment x and the trading strategy tt is defined by 



i=i 



and assumed to satisfy Xq ,7F = x. All the strategies are assumed to be admissible, i.e. 
X^' n > for each t G [0, T] almost surely and self-financing, i.e. 

d 

dXf>* = 4dB t + £ 4dSl t G [0, T]. 

i=l 

A contingent claim, representing future random payoff, is a random variable H > mea- 
surable wrt. Tt- A hedging strategy against H is a pair (x,ir) such that 

P(X T ^ >H) = 1. 

A replicating strategy is a pair (x, tt) such that 

P(X^ ,7r = H) = 1. 

A price of H is defined by 

p(H) := inf {x : 3vr s.t. P(Ay ,7r > F) = l} 
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and, due to the fact that the market is complete, it follows from the general theory that 
p(H) = E[e _rT H], where the expectation is calculated under the measure P. 

If x < p(H) then P(X^' n > H) < 1 for all tt and the question under consideration is 
to find a strategy maximizing the probability of successful hedge, i.e. 

P(X^ >H) — ► max . (2.3) 

We will refer the corresponding function <I>i : [0, +oo) — > [0, 1] given by 

$i(a;) := maxP(4 ,,r > H), 

7T 

as the maximal probability function. The strategy % such that P{X^ > H) = will 
be called the probability maximizing strategy for x. 

We also consider the problem of cost reduction. Let a S [0, 1] be a fixed number 
describing the level of shortfall risk accepted by the trader. Then we are searching for 
a minimal initial cost such that there exists a strategy with the probability of successful 
hedge exceeding 1 — a, i.e. 

x — ► min; 3 vr s.t. P{X^ > H) > 1 - a. (2.4) 

The cost reduction function $2 : [0, 1] — > [0,p(H)] is thus defined by 

$ 2 (a) := min {x : 3ir s.t. P(X^ n > H) > 1 - a) . 

The strategy 7r such that P(X^, 2 ^' W > H) > 1 — a will be called the cost minimizing 
strategy for a. 

In the sequel we study the problem of determining the functions $1 and $2 for the 
contingent claim H of a general form. Then in Section H] specific payoffs are examined. 



3 Main results 

In this section we present a general method of determining functions <3?i and §2- Let us 
start with the auxiliary problems which can be solved via the Neyman-Pearson lemma. 

Assume that we are given two probability measures Pi, P2 with strictly positive density 

and consider two types of optimizing problems 

f Pi\A] — > max, . , 

^ L J 3.5) 
\P 2 [A]<x, 



\ P ^>-'-« (3.6) 
\P 2 [B] — > min, 

where a,x £ [0, 1] are fixed constants. Problem (|3.5p is a classical one appearing in the 
statistical hypotheses testing. Recall, that if there exists a constant c > such that 
P2(jpq>c) = x then the set 

is a solution of (|3.5p . It is not surprising that the solution of the problem (|3.6|) is of a 
similar form. For the convenience of the reader we prove the following. 
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Proposition 3.1 If there exists a constant c > satisfying P\(^p- < c) = 1 — a then the 
set 

^ [dP 2 
\dPi ~ 

is a solution of the problem (|3.6p . 

Proof: Let B be an arbitrary set satisfying P\(B) > 1 — a. We will show that P 2 (B) > 
P 2 (B). The following estimation holds. 



P 2 (B) - P 2 (B) = f(l B - l B )dP 2 = [ (1b - l § )dP 2 + f (1 B - l 6 )m 
>cf (1 B - l S )dPx + c f l B dP x 

= c (J l B dP l - f l§dPi) = < p i( B ) - p i(B)) 



> c[P x (B) - (1- a)) > 0. 



□ 



Let us notice that both optimal sets A, B have a similar form 



l^>c[, (3.7) 



\dP 2 

with suitable constants c > 0. More precisely, for A the constant c is s.t. 

P 2 (g>c)=x (3.8, 

and for B is s.t. 

Pl (S £c ) =1 - a ' (3 ' 9) 

Now, come back to the initial problem of determining functions <J>i, <J>2- Let us start with 
presenting two auxiliary results which are nonrandomized versions of Theorems 2.34 and 
2.42 in 0. 

Theorem 3.2 Let x > 0. If A is a set solving the problem 

{P[A] — > max, 
E[e- rT Hl A }<x, '• U0) 

then <3?i(x) = P(A) and the probability maximizing strategy for x is that one replicating 
the payoff Hl^. 

Let us notice that if x > p(H) then A = Q and thus 3>i(x) = 1. Moreover, if (|3.10|) has a 
solution for every x > 0, then the function $1 is increasing. 
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Theorem 3.3 Let a G [0, 1] be a fixed number. If B is a set solving the problem 

lE[e rl H1 B \ — > mm, 

then $2(0) = E[e _rT Hl^] and the cost minimizing strategy for a is that one replicating 
the payoff HI g . 

Notice that $2(0) = p(H) and if (|3.1ip has a solution for each a E [0, 1] then $2 is 
decreasing. 

Now apply the method of solving the problems (|3.5p and (|3.6p to (|3.1U|) and (|3,lip . 
Notice that f)3.10[) and (|3.1ip can be reformulated to the following form 

( P[A] — > max, 

and 

f P[B] > 1 - a , 
< 1 J - (3.13 
[P*(B) — Mnin, 

where P* is a probability measure given by the density 

dP* H 



dP B[e~ rT H} 

In view of (|3.T[) we are searching for the solutions A, B to (|3.12p . f)3. 13j) in the family of 

sets 

dP ) ( dP dP ) f ~ , c rr l 



dP* ~ J (rfpdP*- j I J ~ E[e-^il] 

where 2y is given by (|2.2p . Denoting, for the sake of simplicity, the constant ^ e - r r H ^ by 
c we see that the optimal sets A, B are of the form 

A c := {Zy 1 > cH} , (3.14) 

where, by (|3.8p and (j3.9fl . c is s.t. 

P*(A C ) = — for A, (3.15) 

and 

P(A C ) = 1 - a for B. (3.16) 
Now define two functions ^1 : [0, +00) — > [0, 1], * 2 : [0, +00) — > [0,p(H)] by 

Vtic) := P{A C ), (3.17) 

* 2 (c) := P*(A C ) ■ E[e~ rT H] = E[e~ rT H1a c ]- (3.18) 

Let us notice that both functions ^1,^2 are decreasing and ^i(0) = 1, ^2(0) = p(H). 
Thus ^2(0) provides the arbitrage free price of the continent claim H. Below we list some 
properties of functions ^1,^2 needed in the sequel. First let us introduce two conditions 
concerning the real function / : R rf — > [0, +00): 
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(CI) X d ({z : f(z) = c}) = for each c> 0, 

(C2) X d ({z : f(z) G (a, b}}) > for each < a < b. 

Above \d stands for the Lebesgue measure on 

Proposition 3.4 a) The function \&i is left continuous with right hand side limits in each 
point of the domain. 



b) The following holds 



lim ^i(c) = P(H = 0). 

c—¥+oo 



Assume that Z T H = f(W T ) where f :R d — > [0,+oo). Then is 

c) continuous if and only if (CI) is satisfied, 

d) strictly decreasing if and only if (C2) is satisfied. 
Proof: a) The function can be written in the form 

* 1 (c) = p(z T H < ^ =F~ ZtH (^-\ c>0, (3.19) 

where F^tH stands for the distribution function of the random variable ZtH . Thus 
has one sided limits for any c > and the left continuity follows from the right continuity 
of Fz t h f° r an y c > 0- Left continuity at c = follows from monotonicity. 

b) The assertion follows from the formula 

*i(c) = P(Z T l >cH\H> 0)P(H > 0) + P(Z T X >cH\H = 0)P(H = 0) 
and the following 

lim PiZZ, 1 > cH I H > 0) = 0. 

c->+oo 1 

c) First show continuity at zero. If c n J, then {Z^ 1 > c n H} n is an increasing family of 
sets and by the continuity of probability we have 



lim fi(c n ) = lim P(Z^ > c n H) = PI I \{Z~ l > c n H}) = PiZ^ 1 > 0) = 1 = ^i(0). 

n 

Taking into account (|3.19|) we see that is continuous for each c > if and only if the 
random variable ZtH = f(Wr) has no positive atoms. In view of the equality 

P{Z T H = c) = P(f(W T ) =c)= C Wt {{z : f(z) = c}) , c> 0, 

and the fact that the distribution of Wt is nondegenerate we see that the continuity of 
\&i is equivalent to (CI). Cy/ T above stands for the distribution of Wt- 
d) For < ci < C2 we have 

-) - P(Z T H <-)= P(f(W T ) G (-, - 
c\J V c 2 / V Vc 2 ci- 



^ 1 {c 1 )-^ 1 (c 2 )=p(ZtH<-)-p(z t H<-^ =p(f(W T )e (-,-]) 
= ^ T ({,:/(,)G(11]}) 



1 1 

- c 2 ' Cl . 

and it follows from the nondegeneracy of the distribution of Wt that the strict monotonic- 
ity of is equivalent to (C 2 ). □ 
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Proposition 3.5 a) The function ^2 is left continuous with right hand side limits in each 
point of the domain. 



b) The following holds 



lim ^(c) = 0. 



Assume that Z T H = f(W T ) where f : R d — ► [0, +00). Then ^ 2 is 

c) continuous if and only if (CI) is satisfied, 

d) strictly decreasing if and only if (C2) is satisfied. 

Proof: a) It follows from monotonicity that one sided limits exist. We show left continuity 
for any c > 0. For c n \ c the family 

{Zj, 1 > c n H} n 

is decreasing and 

Pli^r 1 > c n H} = {H = 0} U {Z^ 1 > cH} = {Z' 1 > cH}. 

n 

Thus by the dominated convergence we have 

lim ^ 2 (c„) = lim n[e- rT Hlr?-i >r J = B[e~ rT H 1 ] = * 2 (c). 

b) For c n t +00 we have 

{Z? 1 > c n H} n I f|{Zf 1 > c n H} = {H = 0} U {Z~ x = +00} = {H = 0}, 

n 

and thus 

J&o* 2 ^ = nl^ [e ~ rTH \^>o n H}\ = ne- rT Hl {H=0} ] = 0. 

c) We show that the right continuity of ^2 is equivalent to (CI). Then continuity follows 
from (a) . For c n \. c > we have 

{Z T X > CnH} | UiZ? 1 > c n H] = {H = 0} U {H > 0, Z T X > cH} 

n 

= {Z- 1 > cH} = {1> cf(W T )}, 

and thus 

^m^Cn) = B[e- rT Hl {1>cf(WT)} }. 

The condition lim n _> +00 * 2 (c) = * 2 (c) holds if and only if P(l > c/(W" T )) = P(l > 
c/(Wt)). The last condition holds for c = and for c > it is equivalent to (CI). 

d) Fix < ci < C2- The inequality 

* 2 (ci) - * 2 (c 2 ) = E[ e - rT Fl { _L </(VKT) <_L } ] > 
holds if and only if P(j- < f(W T ) < jr) > 0. The last condition is equivalent to (C2). □ 
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Now assume that ZtH = /(Wt) for some / : R d — > [0, +00). Let us fix a € [0, 1], 
x > and consider the problem of existence of solutions to the equation 



as well as to 



*i(c) = l-a, (3.20) 



* 2 (c) = x. (3.21) 



In view of Propositions 13.41 and I3"31 it follows that if (CI) is satisfied then ^2 are 
continuous decreasing functions with images (P(H = 0),1] and (0,p(H)] respectively. 
Thus for a G [0,P(# / 0)) and x G (0,p(i?)] the equations (ET20]) and (l3~2TT) do have 
solutions. Moreover, if (C2) is satisfied then the solutions are unique. 

The description of functions <3?i and 3>2 is provided by the following theorem, which is 
the main result of the paper. 

Theorem 3.6 Assume that ZtH = /(Wt) for some f : M. d — > [0, +00) satisfying (CI). 
a) Let c = c(x) G [0, +00) be a solution of the equation 

* 2 (c) = x, x G (0,p(H)). (3.22) 
Then the maximal probability function is given by 



$i(x) 



P(H = 0) for x = 0, 
*l(c(x)) forx£(0,p(H)), 
1 /or x > p(-ff ) 



Moreover, for any x G (0,p(H)) the probability maximizing strategy for x is that one 
replicating the payoff H1a c{x) ■ 

b) Let c = c(a) G [0, +00) be a solution of the equation 

*i(c) = 1 - a, a G [0, P(H / 0)). (3.23) 
Then the cost reduction function is given by 



$2 (a) 



*2(c(o)) /ora€[0,P(#^0)), 
for a G [P{H 7^ 0), 1]. 



Moreover, for any a G [0, P(i? 7^ 0)) i/ie cosi reduction strategy for a is that one 
replicating the payoff H 1a c(q) • 

Proof: The proof is based on the consideration proceeding the formulation of the Theo- 
rem. 

a) If x > p(H) then the hedging strategy is the probability maximizing strategy and 
then clearly $i(x) = 1. Consider the case x G (0,p(H)). By Theorem 13.21 we know that 
&i(x) = P(A), where A is a solution of (|3,10p . The solution of (|3.12p . which is equivalent 
to (pnO"]) . is of the form (I3T41 with c satisfying ([515]) . But (ETT51) is equivalent to <K22\i . 
Thus we have 

$ 1 (x) = P(A c ) = ^i(c), 
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where c is given by the condition ^(c) = x. For x = consider the trivial strategy ir = 0. 
Then P(X^' n > if) = = 0). On the other hand, due to the monotonicity of $i, 

we have $i(0) < lim^ = lim^o *i(c(a:)) = hni 2t+00 * x (z) = P(fT = 0). As a 

consequence we obtain $1(0) = P{H = 0). The second part of the assertion follows from 
Theorem 13.21 

b) If a G 7^ 0),1] then consider a trivial strategy 7r = with zero initial en- 

dowment x = 0. Then X£ ,7r = and thus P{X'^ > H) = P(H = 0) > 1 - a. As a 
consequence we have $2(0) = 0. Now consider the case a G [0,P(P 7^ 0)]. It follows from 
Theorem S3] that $ 2 (a) = E[e~ rT Hl^], where P is a solution to (|3.1ip . The optimal 
solution of (|3.1ip is the same as for (|3.13p and has the form (|3.14p with c satisfying (|3.16p . 
The condition (|3.16p can be written as ^i(c) = 1 — a. Thus we have 

$ 2 (a) = E[e- rT Hl Ac ] = tf 2 (c). 

The second part of the assertion follows from Theorem 13,31 □ 

In virtue of Theorem 13.61 the only problem to determine functions $1, $2 is to find 
functions \l/ 2 and to solve the equations (|3.22j) . (|3.23p . In general, due to the fact that 
^1, ^2 are rather of a sophisticated form, one should not expect to find analytic formulas 
for the constants in (|3.22p . (|3.23|) . However, the equations (|3.22|) . (|3.23p can be solved 
with the use of numerical methods. In the sequel we solve the problem of determining 
functions ^2 for the most common multiple assets derivatives. 



4 Quantile hedging in two dimensional model 

In this section we determine explicit formulas for the functions ^1, ^2 for a few examples 
of popular options. Since our derivatives depend on two underlying assets we simplify at 
the beginning general formulas from Section [3] to the case d = 2. In the calculations we 
base on properties of the multidimensional normal distribution which are recalled in the 
sequel. 

For the case d = 2 we denote the correlation matrix by 



Q 



1 p 
p 1 



Consequently, we have 



Q- 



-1 
p 



p 
-1 



1 + 1 



1 



1 



1 



fl+p 
_L 



+ 



1 



_L 



Hence the density of the martingale measure 



can be written as 



Zj 1 



-A 1 W) r -A 2 W$-BT 



(4.24) 
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where 



. , an - r a 2 - r 
M ■= -3 — 7 + P 



A, 



P — 1 V a l °"2 

ai — r a<i — r 



P — 1 \ CTl CT2 



V 1 + p \/i - fi W 1 + p V 1 - P' a 2 

2 

r ^ 

+ p y/T^pJ a\ ' Vyl + p 1 V 1 - P^ °2 
The formula (|3.14|) for the set A c simplifies to the form 

Ac = \Z T X > CH) = { e MW±+A 2 Wi+BT > cH j ; 

and consequently formulas (|3.17p . f|3. 18|) become 

tt 2 (c) = E[ e - rT J Efl {eAlW i +A2 ^ +ST ^ cH} ]. 



Now set the notation concerning the multidimensional normal distribution and recall 
its basic properties, which can be found in standard textbooks on probability theory or 
statistics, see for instance [8] . A random vector X taking values in M. d has a multidimen- 
sional normal distribution if its density is of the form 

f x {x) = r • e -|(^-™) T s- 1 (^-™) ) x e R d (4.25) 

(2vr)2(detS)2 

where m 6 R rf is a mean of X and X is a symmetric positive definite d x d covariance 
matrix of X. The fact that X has a density ()4.25p will be denoted by X ~ Nd(m, S) or 
£(X) = Nd(m, £). If <i = 1 then the subscript is omitted and N(m, a) denotes the normal 
distribution with mean m and variance a. If X ~ N^(m, S) and A is a x c! matrix then, 

AX ~ N k (Am, AT,A T ); (4.26) 

in particular if a 6 M rf then 

a T X ~ iV(a T m, a T Sa) . (4.27) 

Let X be a random vector taking values in R d and fix an integer < k < d. Let us divide 
X into two vectors X^ and X^ with lengths k, d — k respectively, i.e. 

X^ = (Xi, X 2 , Xk) T 1 X^ = (Xk+i, Xk+2, Xd) T ■ 

Analogously, divide the mean vector m and the covariance matrix S 



m ( 2 ) 



S (21) S (22) 
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sothatEXW = m«, EX( 2 ) = , CovX^ = E( n ), CovX^ = E( 22 \ Cov{XW, X^) = 
S (12) = s (2i) T Denote by £ (^W | X^ 2 ) = x^) the conditional distribution of X^ given 
jsf(2) = x (2) £ If s (22) ig noris i n g U i ar then 

£ (x« | X^ = x ( 2 )) = N k (m^(x^)^ u \x^)), (4.28) 

where 

m«(x( 2 )) = m« + S( 12 )s( 22 ) _1 (x( 2 ) - m( 2 )), 

E^ 1 )^ 2 )) = E< u > - S( 12 )S( 22 ) _1 E( 21 ). (4.29) 

Actually the conditional variance E^ 11 ) (x^ 2 ^) does not depend on x^ but we keep the nota- 
tion for the sake of consistency. The conditional density will be denoted by /x( 1 )|x( 2 )=a;( 2 ) 
where x^ G In particular if (X, Y) is a two dimensional normal vector with parame- 
ters 

OH <7l2 



mi 
m 2 



C21 022 



then 



where 



£(X I Y = y) = N(m 1 (y),a 1 (y)), 



(J C7 2 

m 1 (y) := m 1 -\ (y - m 2 ), o"i(y):=crn — . (4.30) 

C22 CT22 

If X is a random vector then its distribution wrt. the measure P will be denoted by C(X) 
and its density by fx- Analogously, /x(i)ix( 2 )=a:( 2 ) ( x ) stands for the conditional density 
with respect to the measure P. 

In the following subsections we will use the universal constants: Ai,A2,B defined in 
(|4.24p as well as oj, a 2 , b, 5i, 02, b introduced below. 
Fix a number K > 0. One can check the following 

{S T > K } = {Wt > 01} = {w T > «i} > ( 4 - 31 ) 

{Sf > if} = {W| > a 2 } = > 5 2 } , (4.32) 

{S T > S T } = {a x W T - a 2 W$ > b} = {<j x W t - a 2 W% > £>} , (4.33) 

where 

ai := ^ (in | - («i - ~^)t) , 3l := i (in * - (r - ia 2 )T 
a 2 := ^ (in * - (a 2 - ^2)t) , a 2 := 1 (in J - (r - 

6 := In (|) + (a 2 - Ox - ^1 - a 2 ))T, b := In (|) - \{a\ - a\)T. 

In all the formulas appearing in the sequel it is understood that In = —00 and $ stands 
for the distribution function of N(0, 1). 
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4.1 Digital option 

In this section we determine V&i, ^f 2 for the payoff 

L {Si>S'2}> 



(4.34) 



By (jUgSD we have 

*!(c) = P(A C ) = P (^+^w?+«r > cPl {s i > s , ^ 

= P (AiWj + v4 2 Wf + PP > ln(cP), 4 > S 2 T ) + p( e AiW± + A 2 w*+BT >q iS i < s ^ 

= P (A-lW^ + A 2 W% + BT> \tl(cK) I C7i - a 2 W£ > b) P{a x W\ - a 2 W} > b) 

+ P{a x W^ - a 2 W% < b). (4.35) 

Let us notice that 

X := 



A X W\ + A 2 W% ' 




'At A 2 ~ 






o\W\ - a 2 W^ 




o\ -cr 2 







so in view of (|4.26p we have X ~ A^(0, S), where 

(Ai + A 2 r)TAi + (Air + A 2 )TA 2 (ax - a 2 r)TA 1 + (a x r - a 2 )TA 2 
(a \ - o 2 r)TA\ + (a\r - a 2 )TA 2 (a\ - o 2 r)To\ - (p\r ~ o 2 )To 2 
In virtue of (|4.30[) we have 

L(A X W\ + A 2 Wl I a x W\ - o 2 W% = y) = N(m(y),a(y)), 

where 



m(y) = y 



(a 1 - a 2 r)Ai + (a\r - a 2 )A 2 _ 
(<7i - o 2 r)o\ - (o\r - a 2 )a 2 ' 



T(A 1 a 2 + A 2 a 1 ) i (p l -l) 
—a\ + 2po\o 2 — a\ 



By gJ7D we have: 01 W4, - ct 2 W| ~ iV(0, T{a\ - 2pa 1 cr 2 + erf)). Going back to (g^5]) we 
have 

/■+00 

*i(c) - / P(AiWf + A 2 W| > ln(cP) - PP I a x W^ - <t 2 W* = y)L lW ^-a 2 w-(y)dy 

J b 

+ P(a x W\ ~ cr 2 Wf < b) 



" » / + PP - ln(cP) \ 

V v^) ) f ^-^ {v)dy + ^ 
Now let us determine \E , 2- In virtue of (j4.33p we have 
* 2 (c) = e - rT E[Pl Ac ] = e" rT E[Pl 



A /T((7if-2^7iCT 2 +<^) y 

(s£ > Sl,Z T l > 1 K \ 



e~ rT KP 



{ S rj-i 5* rji } 



-rT 7 



#P ( Z^ 1 >cK\S^> Sf ) P{S^ > s% 



e-* r KP(e AlW ' +A * w * +BT > cK \ a,W^ - a 2 W* > fyPfaW} - a 2 W* > b) 

A+oo 



+00 



K ■ / P(A{W^ + A 2 W% > ln(cP) + A 1 — — -T + A 2 — — -T - PP | a x W\ - cr 2 W| = y) 



ot 2 - r r 



0"! 



(To 



(y)dy 



e~ rT K 



+~ ^ m{y)-HcK)-A l ^T-A 2 ^T + BT s 



14 



o- 1 w r i,-<T 2 H'; 



4.2 Quantos 

4.2.1 Quanto domestic 

The contingent claim is of the form 

H = S^(S^ - K) + , K>0. (4.36) 
At the beginning let us notice that 

^wI+AiW'+bt > cS 2(gi _ K) j = | (A2 _ a2)W 2 > „( CjW ^)} = |(A 2 - a 2 )WZ > w(cM)} , 

(4.37) 



where 



(c,x) := In (^e^-^'-B^-Aix^ig^-ia^T+^x _ ^ 



u>(c, x) := In 



By (OTD and (14371) we have 

+ P LAiW}+A a W}+BT > | 5^ < if) P(4 < tf) 

= P ((A 2 - o- 2 )Wf > In (c^e (Q2 -3-2-s)T-A 1 ^ (5 i _ ^ i w i > ffl ^ p ^ w i > ai j + < ai) 

= ^ + °° P((A 2 - a 2 )W| > u(c, W£) | W£ = + * 

The conditional distribution is given by 

C((A 2 - a 2 )W% | W£ = a?) ~ 7V(m(x),<j(x)), 
where m(a;), <r(a;) are given by (|4.29p . Hence we have 



*i(c) 



To avoid technicalities assume that A 2 ^ o 2 . We have 



* 2 (c) = e rT E |S T (S T - ^) + l {e A lW i+4 2l f| + BT >cS 2 (s ,_ 



K) + } 



e- rT E 



S T (S T K ) + 1 {e A 1 W^ + A 2 W^, + BT^ cS ^ {slT _ K)+} | ■';/ 



< K 



(5i-if)+} 



S£ > if 



S T {S T K)l AlW i +A2W 2 +BT 



>cSUSi.-K)} \ S T>K 



P(S^ < K) 
P{S^ > K) 
P{S^ > K). 
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By (jOTD and (Q71> we have 

* a (c) = e~ rT E [flgefr-W^^^Cr-W^^ _ K)!^^^^ \ W\ > 3i] > Sx) 



-rT 



e " / E 



{( J 4 2 - CT2 )H'|>«;(c,Wi)} I W4 - X 



^ e (r-**I)T+<r2V^/£l e (r-£a?)T+ffiW£ _ 



(c,x) 



+ 00 /"+00 



/•-t-oo /> 

^J^V 1 -^ 7 / / 

«/ di ■/ tt 



■t-— ^ — y 

A 2~ cr 2 



Si */ u;(c,a;) 



f { A 2 -a 2 )w^m= x (y) d yfw^ x ) dx - 



It follows from that £((A 2 -<7 2 )Wf | = x) = N((A 2 -a 2 )px, T(l- p 2 )(A 2 -a 2 ) 2 ) 

and hence 



* 2 (c) = e- rT SlSle^ r ---<-^i 



- e- rT KS^ r -i^ T 



^2nT{l-p^){A 2 -a 2 f 



I {y-(A<2-<T 2 )pxr 



^27^(1 -^(Aa-tra) 2 Ax A (c ,x) 



00 />+oo 



(c,a:) 



CT 2 ^ 1 (^-{^2~ ff 2)P^) 



4.2.2 Quanto foreign 

The payoff is of the form 



H=[Sk 



'T 



K 

Srp 



, K > 0. 



First let us notice that 
K 



(4.38) 



where 



d:=ln 



Ql o2 

DqDq 



ai + a 2 --(cr 1 +cr 2 ) T, rf:=]n-j- 



01 o 



2 v 2r--(^ + a 2 2 )]r, 



and 



,A 1 W^+A 2 W^+BT > c ( 5I 



7 r 



{A-lWt + (A 2 + a 2 )W% > v(c, a x W\ + a 2 W%)} 

{a x W^ + (A 2 + o- 2 )Wf > u/(c, trjWj + a 2 W$)} , 

(4.39) 



where 



W ( C)2 ) : =ln^ e (^- Q2 - B ) T (5 1 5 V 1+Q2 -5(-?+-2 2 )^_^^ , 

B (C*) :=ln [j-e^^^* 2 ^ (Sl%e^\rt+®F+* - K 



16 



By (OHJ) we have 



*i(c) - P ( e ^^+A 2 ^+BT > c ( <?1 _ |L 

iJrp 



= P ( eAiW*+A a W|+BT > c l S l_K\\ aiW l + U2W 2 > j \ p ((Ji W l + U2W 2 > rf) 



Srp 



+ P(a 1 W^ + a 2 W^ < d). 
As a consequence of (|4,39p we obtain 

*i(c) = P(AiWf + (A 2 + (7 2 )W| > u(c, ctiW£ + <t 2 W|) | a x W\ + o- 2 W| > d)P{a x W\ + <t 2 Wt 

+ P((TiWj + <7 2 W| < d) 
/•-foo 

= / P(AiW4 + (A 2 + a 2 )W| > i;(c, z) | <nW£ + a 2 W| = z)f aiW i +a2W 2{z)dz 

J d 

+ P(criW4 + °2 W T < d )- 
By (jllSUD we have 

L{A X W\ + (A 2 + o- 2 )W| | oiWr + o- 2 W"| = *) = N{m{z), <r(z)), 



where 



(Ai + (A 2 + a 2 )p)a 1 + (Aip + A 2 + a 2 )c7 2 

mlz) := o o > 

erf + 2pa x a 2 + <t| 

ff(z) := T j (Ai + (A 2 + a 2 ) p) A l + (A lP + (A 2 + a 2 )) (A 2 + a 2 ) 

{{Ax + (Ag + a 2 ) p) oi + (Aip + (Ag + <r 2 )) <g) " 
of + 2po-icJ 2 + of 



and thus 
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By (PHj) and (1P9D we have 



* 2 (c) = e~ rT E 



e- rT E 



S T - -g2 j 1 {e ^l"'i + ^2"'2 + f5T >c(s ^__^_ )} I CTlW4 + 0-2 W| > d 



P{<T X W?t + a 2 W% > d) 



Sn ■ 



K 

J ^{A 1 W^ + (A 2 +a 2 )W 2 >w(c,a 1 W^+a 2 W 2 )} 



(JiWt + (T 2 Wt > d 



P(a- 1 W^ + a 2 WT > d) 



e~ rT I E 

Id 



S T 1 {A 1 W± + (A 2 +a 2 )W2>w(c,<7 1 W±+a 2 W2)} I °lW T + <T 2 W T - Z f aiW i +rT2W 2(z)dz 



r+00 

e~ rT K / E 

Jd 



'g2 1 {A 1 W^ + (A 2 +a 2 )W^>w{c^ 1 W^+cr 2 W 2 )} I a ^ W T + a * W T ~ z 



Using ()4,28p we find the conditional distributions 

£(Wt,AiW^ + (A 2 + a 2 )WT I (?iW^ + a 2 W% = z) = N 2 (M 1 (z), S 1 ^)), 
C(W$, AiW^ + (A 2 + <7 2 )Wr I aiV^f + ct 2 ^| = z) = N 2 (M 2 (z), S 2 (z)), 

where M 1 (2;), M 2 (z), S 2 (z) are determined by (|4.29p . As a consequence we obtain 



p+00 p+00 p+00 
* 2 (c) = e-^S^e^-^?) 7, / / / e^ x F\x,y)dy dxf] 

Jd J — 00 J w(c,z) 



[z)dz 



-r T £ e -{r-\a\)T 
°0 



OO /"+00 /" + OO 



00 Jw(c,z) 



e ^ x F 2 (x,y)dy dxf aiW , +(72W ,(z)dz, 



where F 1 , F 2 stand for the density functions of the two dimensional normal distributions 
N 2 {M 1 (z),E 1 (z)), N 2 {M 2 (z),Y?{z)) respectively 

4.3 Outperformance option 

The problem is studied for 

H = (max{5£, Sf.} - K) + , K > 0. 

Let us notice that 

{ e AiWi+A 2 W*+BT > c(5 l _ j = | AaW ,2 > Wl ( Cj ^)} = |^ 2 M?2 > ^ (440) 

{e A ^ +A > w * +BT > c(S* - K)} = {A 1 W^>v 2 (c,W})} = [a x W^ > w 2 (c, w|)} , (4.41) 
where 

/ 1 9> m . "1 r ~ / 1 9s m . 

A 2 y - BT, 



«l(c,x):=ln C^e^-S^i)^!^ _ ^ _ fiT U2 ( C)2/ ). = l n CiS r2 e (a 2 - 5 <r 2 )T+ CT3!/ 

Wl (c,x) := In (ce-^xx+^x^F-flJT^igCr-^jr+a.x _ x) ) 



012 — T „ 

+ — T, 



w 2 



(c,y) := In (t 



,-A 2 y+(A 2 !^~B)T^ s 2 e (r-±a 2 2 )T+ a2 y _ 



K) 



u 2 

ai — r 



T. 
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By (|OTj) . (Q2I) . (Q3|) we have 

*!( C ) = P ^iWj+^Wj+BT > c(5 l v ^2 _ K) + ^ 

=P f e AiWi,+A 2 W*+BT > _ jjQ | ^1 > au(JlW l _ a2W 2 > ^ p^ > ^^^1 _ > ft j 

+ P ^Wt+^b't+bt > I Wf < ai,criW4 - 0-2 W| > b) P(W% < a u a{W^ - a 2 W% > b) 

+ P > o I Wt < o 2 ,criWf - cr 2 W| < P(W| < a 2 ,aiW£ - cr 2 W| < 6). 

By IfOQjL (gHH we have 

*i(c) = P (A 2 Wr > «i(c, W£) | > ai,criWj - (7 2 W| > 6) P(W£ > at,a{W^ - cr 2 W| > b) 
+ P(M4 < ai, (tiWt - cr 2 W| > 6) 

+ P (ii^ > u 2 (c, Wf ) | W| > a 2 , o-i W4 - cr 2 W| < 6) P(W| > a 2 , cti w£ ~ °2W| < 6) 

+ P(W| < a 2 , (71 - C7 2 W| < 6), 

Let mi(y, z),m 2 (x, z),a±(y, z), a 2 (x, z) be the means and variances of the conditional dis- 
tributions 

C(AxW£ | WZ = y,a 1 W£-a 2 WZ = z)=N(m 1 {y,z),a 1 (y,z)), 

£(A 2 W% | W^ = x,a 1 W±-a 2 W% = z)=N(m 2 (x,z),a 2 (x,z)), 
given by (|4.29p . Then we have 

= / + °° f W ^ iW ,_ a2W ,{ X , Z )dzd X 

Jai Jb \ y/(T 2 {x,z) J 

ai z'+oc 

fw^,a 1 w^-<7 2 wz ( x , z)dzdx 



— GO 



+ °° /" 6 $ | m^y.z) - v 2 (c,y) 
a 2 J-oo \ yjai{y,z) 



fw^,,a 1 Wi,-<T 2 W?,(y^ z)dzdy 



a 2 fb 

fwi^w^-aiw^iyi^dzdy. 



oo */ -co 



By (|OI]> . (|Q2]> . flES}, P0| . (IOTP we have 
* 2 (c) = e--E ((# V 4 - if) + l {e ^ + ^^ c(m+}/ 

= e~ rT E ((S^ - ^)l {e A 1 »'i + A 2Bf | + BT >c(s ^ Jf)} I St > K, S T > St J P(St > K, S T > S T ) 

+ e~ rT E ((Sf - if)l {eAlW i +A2W 2 +BT>c( _^ )} | Sf > K, S T < SfJ P(Sf > K, S T < Sf) 

= e- rT E ((S* T - K)l {A2W ^ wiic W i )} | W\ > Sx.trxW^ ^ <7 2 W? > ft) P(W T > a^W^ - a- 2 Wf > 6) 

+ e- rT E - if)l {Al f>> W2(Ci ^ )} I W| > a 2 , (7i - (x 2 ^ < ^ p(w?2 > _ CT ^2 < ^_ 
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Let mi(y, z),a±(y, z) and m 2 (x, z), a 2 (x, z) denote means and variances of the conditional 
distributions 

t[A x W\ I Wla 1 W^-a 2 W^j = N ( mi (y, z), a^y, z)) , 

c(a 2 W% I W±,a 1 W£-a 2 w£) = N (m 2 (x, z),a 2 (x, z)) . 
Finally we obtain 

- c r ^'- k "' )T+ '" - «) * ( m,( ^ i> ) 

4.4 Spread option 

The payoff is of the form 

H= (S^-S^ -K) + , K>0. 

For any y G R 

- S^-t^+^j/ _ ^ > 0} = {H4 > e(y)}, (4.42) 

where 

-(--rO T )' 

By conditioning wrt. we obtain 

/+oo 
p (>^+A 2 ^ +S T > C(5 1 _ <.2 _ x)+ ]W 2 =y \ fw , [y)dy 
-oo 

/+oo . . 

p (>^ + A 2 ^ +B T > C(5 1 _ <.2 _ ^1 > e(y) ! ^2 = \ 
-oo 

/'+00 . , 

p U^+^WJ+BT > Qj W l < e{y) \ W 2 =y \ fw , {y)dy 
-oo 

/+oo 
-co 

/+oo 
P(w£ <e(y)| W| = y)/ W r,(»)di/. 
-oo 

Let 

£(W±\W* = y) = N(m(y),a(y)). 
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e(y) := — In 

0"! 



i_(52 e («2-ia|)T + a 22/) 



Then we have 

= / + °° / f wilw > =y (x)dxf w *(y)dy + / + °° $ f ffc^ ] /^( y )dy, 

J-oo J5( C ,y)n(e(j/),+oo) ./-oo \ V "^) / 

where 5(0, y) is a set defined y G R by 

5(c,y) := {x : ^v+bt^ _ ^i^-i^T^x > _ c ^ e ( aa -^yr+ ffav + K) | 

For the practical applications it is necessary to find a closed form of the set S(c,y). In 
the formulation of the next result we will use the solutions of the equation 

g(x) = -c(S%e( a2 -^ T+a2V + K), (4.43) 

where g{x) := e A2V+BT e A±x —cSQe^ ai ~^ a ^ T e aiX . These solutions can be found numerically. 

Proposition 4.1 The set S(c,y) is of the form 
a) if A\ > cf\ and 

(i) ifg(x) > -c(S^ a2 -^ T+<72y + K) then S(c,y) = (-oo,+oo), 

.1^ 



(ii) if g{x) > -ciS^e^-^T+^y + K ) then s ^ y ^ = (_oo,xi) U (x 2 ,+oo), where 
x\ < X2 are the unique solutions of (|4.43p . 



Above, x stands for — ^h- In ( — Me 2V+ \ _ 
b) if A\ = o\ and 



(i) e M y +BT > c5 1 e( ai -5 CT i) T then S(c,y) = (-oo,+oo), 

(ii) e A 2y+BT cSQe( ai ~2 a ^ T then S(c,y) = (— oo, xq), where xq is a unique solution 

of 

c) if A\ < o\ then S(c,y) = (— oo,xo), where xq is a unique solution of (|4.43p 

Proof: a) One can check that g has a minimum at the point x and is decreasing on 
(— oo,x) and increasing on (x,+oo). Hence (i) and (ii) follow. 

b) The formulas for S(c, y) follows from the simplified form of the function g(x) = 

( e A2 V +BT _ c Sl e K-H)T) e A ia; _ 

c) It can be checked that g is strictly increasing on the set {x : g(x) < 0} and that 
lim. r _ 5 . +00 g(x) = — oo. Thus (|4.43p has a unique solution and the form of the set S(c,y) 
follows. □ 

Now let us determine $2- One can check that for y G R 

{S^ - S$e (r -^ )T+<72y - K > 0} = {W£ > f(y)}, (4.44) 

where 

^ p (r-i CT 2)T+<T 2 j/N 



f(y) := - ( In 
01 



For y G R define 
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We have 

* a (c) = e- rT E [(S T - S T - K)+l {eAlW , +A2W , +BT _^ )+} 



/+oo _ 
E (S T -S t - K) + l AlW i +A2W 2 +BT> c(s1 _ S 2_ k)+} \W T = y f^2(y)dy 

/+oo _ ^ 
E [(fii -S T - ^\ 6 ^ + ^ + BT> <a i_^_^ } \^> m} \W T = y\ f w ,(y)dy 



e~ rT I E 



+00 

(s T ~ s T - k) + i AlW i +A2W 2 +BT , 1 w T = y fw*{y)dy 



e- rT 



+00 

■00 JS(c,y)n(f(y),+oo) 



{e A 1 ^ + A 2 ^ + B T> c(s ^_ s 2_ A - ) + } J - {V l/^</(y)} 

^(r-^JT+o!- _ figefr-W^y _ ^ f~, m ^( x ) dx f~,(y)dy. 



The explicite form of the set £(c, y) can be established in the same way as for S(c, y) in 
Proposition 14.11 
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